Light clusters production as a probe to nuclear symmetry energy 
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Using an isospin-dependent transport model for heavy-ion collisions induced by neutron-rich nuclei 
at intermediate energies, we study the production of light clusters such as deuteron, triton, and 3 He 
via coalescence of nucleons. We find that both the yields and energy spectra of these light clusters are 
affected significantly by the density dependence of nuclear symmetry energy, with a stiffer symmetry 
energy giving a larger yield. 
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The energy per particle of an asymmetric nuclear 
matter with density p and an isospin asymmetry 5 — 
(p n — p p )/p, where p n and p p are, respectively, its neu- 
tron and proton densities, is usually approximated by a 
parabolic law Q, i.e., 



E(p,6) = E(p,0) + E sym (p)S 2 , 



(1) 



where E{p, 0) is the energy per particle of symmetric nu- 
clear matter and E sym (p) is the nuclear symmetry energy. 
Although the nuclear symmetry energy at normal nuclear 
matter density po = 0.16 fm -3 has been determined to 
be around 30 MeV from the empirical liquid-drop mass 
formula 0, its values at other densities are poorly 
known. Studies based on various theoretical models also 
give widely different predictions % . Lack of this knowl- 
edge has hampered our understanding of both the struc- 
ture of radioactive nuclei [H El, 13, IH and many important 
issues in nuclear astrophysics [|J lift lllj | , such as the nu- 
cleosynthesis during pre-supernova evolution of massive 
stars and the properties of neutron stars |H.lll). However, 
recent advance in radioactive nuclear beam facilities pro- 
vides the opportunity to study the density dependence of 
the nuclear symmetry energy. Theoretical studies have 
already shown that in heavy-ion collisions induced by 
neutron-rich nuclei, the effect of nuclear symmetry energy 
can be studied via the pre-equilibrium neutron /proton 
ratio 01 > the isospin fractionati on Il3l H3. HH Hq , the 
isoscaiing in multifragmentation [17| . the proton differ- 
ential elliptic flow [18J, the neutron-proton differential 
transverse flow 0], the tt~ to tt + ratio , and two- 
nucleon correlation functions [2lJ . 

In this work, we study the production of deuteron, tri- 
ton, and 3 He in heavy-ion collisions induced by neutron- 
rich nuclei by means of the coalescence model based on 
the nucleon phase space distribution functions from an 
isospin-dependent transport model. It is found that both 
the multiplicities and energy spectra of these light clus- 
ters are sensitive to the density dependence of nuclear 
symmetry energy but not to the isospin-independent part 
of nuclear equation of state and the in-medium nucleon- 
nucleon cross sections. Therefore, light clusters produc- 



tion in heavy-ion collisions induced by neutron-rich nuclei 
provides another possible method for extracting useful 
information about the nuclear symmetry energy. 

The coalescence model has been used extensively in 
describing the production of light clusters in heavy-ion 
collisions at both intermediate [23, |23j, |2J, |25j and high 
energies [jfl Et) ] . In this model, the probability for pro- 
ducing a cluster is determined by its Wigner phase-space 
density and the nucleon phase-space distribution at freeze 
out. Explicitly, the multiplicity of a M- nucleon cluster 
in a heavy ion collision is given by |26l | 



Nm = G J dr il dq il ---dr iM _ 1 dq iM _ 1 

x ( P7( v ii,<lii •••r iM _ 1 ,q iM _ 1 )).(2) 

il>i2>--->iM 

In the above, r ix , • • • , r iM _ 1 and q ix , • • • , qi M _ t are, re- 
spectively, the M — 1 relative coordinates and momenta 
taken at equal time in the M-nucleon rest frame; p^ is 
the Wigner phase-space density of the M-nucleon cluster; 
and (• • • ) denotes event averaging. The spin- isospin sta- 
tistical factor for the cluster is given by G, and its value is 
3/8 for deuteron and 1/3 for triton or 3 He, with the lat- 
ter including the possibility of coalescence of a deuteron 
with another nucleon to form a triton or 3 He 28]. 

For the deuteron Wigner function, it is obtained from 
the Hulthen wave function, i.e., 



a/3(a + (3) e" 
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2n(a - (if 



(3) 



with parameters a = 0.23 fm -1 and (3 = 1.61 fm -1 to 
reproduce the measured deuteron root-mean-square ra- 
dius of 1.96 fm [2^|. As in Ref. (2(|, we express the 
Hulthen wave function in terms of 15 arbitrary Gaus- 
sian functions and determine their strengths and ranges 
by a least square fit. For triton and 3 He, their Wigner 
functions are obtained from the product of three ground- 
state wave functions of a spherical harmonic oscillator 
with its parameters adjusted to reproduce the measured 
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root-mean-square radii of triton and 3 He, i.e., 1.61 fm 
and 1.74 fm, respectively [3(J. Normal Jacobian coordi- 
nates for a three-particle system are then introduced to 
derive the Wigner functions for triton and 3 He as in Ref. 

m 

The space-time distribution of nucleons at freeze out is 
obtained from an isospin-dependent Boltzmann-Uehling- 
Uhlenbeck (IBUU) transport model (e.g., [H El El 
For a review of the IBUU model, we refer the 
reader to Ref. [l|. The IBUU model includes explicitly 
the isospin degree of freedom through different proton 
and neutron initial distributions as well as their different 
mean-field potentials and two-body collisions in subse- 
quent dynamic evolutions. For the nucleon-nucleon cross 
sections, we use as default the experimental values in free 
space. In order to study the effects due to the isospin- 
dependence of in-medium nucleon-nucleon cross sections 
^medium) we also use a parameterization obtained from 
the Dirac-Brucckner approach based on the Bonn A po- 
tential 32]. For experimental free-space cross sections, 
the neutron-proton cross section is about a factor of 3 
larger than the neutron-neutron or proton-proton cross 
sections. On the other hand, the in-medium nucleon- 
nucleon cross sections used here have smaller magnitudes 
and weaker isospin dependence than er cxp but strong den- 
sity dependence. For the isoscalar potential, we use as 
default the Skyrme potential with an incompressibility 
A'o = 380 MeV. This potential has been shown to re- 
produce the transverse flow data from heavy ion colli- 
sions equally well as a momentum-dependent soft poten- 
tial with K = 210 MeV [HQ. 

The IBUU model is solved with the test particle 
method j3{| . Although the mean-field potential is evalu- 
ated with test particles, only collisions among nucleons in 
each event are allowed. Light cluster production from co- 
alescence of nucleons is treated similar as nucleon-nucleon 
collisions, i.e., only nucleons in the same event are al- 
lowed to coalesce to light clusters. Results presented in 
the following are obtained with 20, 000 events using 100 
test particles for a physical nuclcon. 

For the density dependence of sym metry energy, we 
adopt the parameterization used in Ref. [3j| for studying 
the properties of neutron stars, i.e., 



E sym (p) = E sym (p ) ■ u 1 , 



(4) 



where u = p/po is the reduced density and E sym (po) = 30 
MeV is the symmetry energy at normal nuclear matter 
density. In the following, we consider the two cases of 
7 = 0.5 (soft) and 2 (stiff) to explore the larg e range of 
E sym (p) predicted by many-body theories |10j . 

We consider the reaction of 52 Ca + 48 Ca, which has 
an isospin asymmetry S = 0.2 and can be studied at fu- 
ture Rare Isotope Accelerator Facility. In the present 
study, nucleons are considered as being frozen out when 
their local densities are less than po/8 and subsequent 
interactions do not cause their recapture into regions of 
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FIG. 1: Time evolutions of the average multiplicity of (a) 
deuteron, (b) triton, and (c) 3 He from central collisions of 
52 Ca+ 48 Ca at E = 80 MeV/nucleon by using the soft (solid 
curves) and stiff (dashed curves) symmetry energies with a 
stiff nuclear compressibility A'o = 380 MeV and free nucleon- 
nucleon cross sections. Results using the soft symmetry en- 
ergy and free nucleon-nucleon cross sections but Kq = 201 
MeV are shown by dotted curves, while those from the soft 
symmetry energy and Kq = 380 MeV but in-medium nucleon- 
nucleon cross sections are given by dash-dotted curves. 



higher density. Shown in Figs. H ( a ), (b) and (c) are 
time evolutions of the average multiplicity of deuteron, 
triton, and 3 He from central collisions of 52 Ca + 48 Ca at 
E = 80 MeV/nucleon by using the soft (solid curve) and 
stiff (dashed curve) symmetry energies. It is seen that 
production of these light clusters from a neutron-rich 
reaction system is very sensitive to the density depen- 
dence of nuclear symmetry energy. Final multiplicities of 
deuteron, triton, and 3 He for the stiff symmetry energy 
is larger than those for the soft symmetry energy by 51%, 
73%, and 100%, respectively. This is due to the fact that 
the stiff symmetry energy induces a stronger pressure in 
the reaction system and thus causes an earlier emission of 
neutrons and protons than in the case of the soft symme- 
try energy [3j| , leading to a stronger correlations among 
nucleons. Furthermore, the soft symmetry energy, which 
gives a more repulsive symmetry potential for neutrons 
and more attractive one for protons in low density region 
(< Po) than those from the stiff symmetry energy, gen- 
erates a larger phase-space separation between neutrons 
and protons at freeze out, and thus a weaker correlations 
among nucleons. The larger sensitivity of the multiplic- 
ity of 3 He to the nuclear symmetry energy than those of 
triton and deuteron as seen in Fig. ^ reflects the fact 
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that the symmetry energy effect is stronger on lower mo- 
mentum protons than neutrons IB?! . 

Although the density dependence of nuclear symme- 
try energy affects appreciably the yield of light clusters, 
changing the incompressibility from Ko = 380 to 201 
MeV (dotted curves) or using cr me dium instead of er oxp 
(dash-dotted curves) only leads to a small change in the 
yield of these clusters as shown in Fig. ^ This implies 
that the relative space-time structure of neutrons and 
protons at freeze out is not sensitive to the equation of 
state (EOS) of symmetric nuclear matter and in-medium 
nucleon-nucleon cross sections. 
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FIG. 2: Kinetic energy spectra in the center-of-mass system 
for (a) deuteron, (b) triton, and (c) 3 He from central collisions 
of 52 Ca+ 48 Ca at E = 80 MeV/nucleon by using the soft (solid 
squares) and stiff (open squares) symmetry energies with a 
stiff nuclear compressibility Ko = 380 MeV and free nucleon- 
nucleon cross sections. 

The kinetic energy spectra in the center-of-mass sys- 
tem for deuteron, triton, and 3 He are shown in Fig. 
for both the soft (solid squares) and stiff (open squares) 
symmetry energies. It is seen that the symmetry energy 
has a stronger effect on the yield of low energy clusters 
than that of high energy ones. For example, the sym- 
metry energy effect on the yield of deuteron, triton, and 
3 He is about 60%, 100%, and 120%, respectively, if their 
kinetic energies are around 10 MeV, but is about 30%, 
40%, and 85%, respectively, if their kinetic energies are 
around 100 MeV. This follows from the fact that lower 
energy clusters arc emitted later in time when the size 
of nucleon emission source is relatively independent of 
nuclear symmetry energy, leading thus to a similar prob- 
ability for nucleons to form light clusters. Since there are 
more low energy nucleons for a stiffer symmetry energy, 



more light clusters are thus produced. On the other hand, 
higher energy nucleons are emitted earlier when the size 
of emission source is more sensitive to the symmetry en- 
ergy, with a smaller size for a stiffer symmetry energy. 
The probability for light cluster formation is thus larger 
for a stiffer symmetry energy. This effect is, however, re- 
duced by the smaller number of high energy nucleons if 
the symmetry energy is stiffer. As a result, production of 
high energy light clusters is less sensitive to the stiffness 
of symmetry energy. This is different from that seen in 
the correlation functions between two nucleons with low 
relative momentum, where the symmetry energy effect 
is larger for nucleon pairs with higher kinetic energies 
[2ll | , as they are only affected by the size of the emission 
source, not the number of emitted nucleon pairs. 




FIG. 3: The ratio t/ 3 He as a function of the cluster kinetic 
energy in the center-of-mass system by using the soft and stiff 
symmetry energies. The lines are drawn to guide the eye. 

The isobaric yield ratio t/ 3 He is less model-dependent 
and also less affected by other effects, such as the feed- 
back from heavy fragment evaporation and the feed-down 
from produced excited triton and 3 He states. In Fig. [31 
we show the t/ 3 He ratio with statistical errors as a func- 
tion of cluster kinetic energy in the center-of-mass system 
for the soft (solid squares) and stiff (open squares) sym- 
metry energies. It is seen that the ratio t/ 3 He obtained 
with different symmetry energies exhibits very different 
energy dependence. While the t/ 3 He ratio increases with 
kinetic energy for the soft symmetry energy, it decreases 
with kinetic energy for the stiff symmetry energy. The 
symmetry energy thus affects more strongly the ratio of 
high kinetic energy triton and 3 He. For both soft and 
stiff symmetry energies, the ratio t/ 3 He is larger than the 
neutron to proton ratio of the whole reaction system, i.e., 
N/Z= 1.5. This is in agreement with results from both 
experiments and statistical model simulations for other 
reaction systems and incident energies [38l |3^, 0, |4l| . 
It is interesting to note that although the yield of lower 
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energy triton and 3 He is more sensitive to symmetry en- 
ergy than higher energy ones, as shown in Fig. [21 their 
ratio at higher energy is affected more by the symmetry 
energy. Moreover, the energy dependence of the ratio 
t/ 3 He is insensitive to the EOS of symmetric nuclear mat- 
ter and in- medium nucleon-nucleon cross sections. These 
features thus imply that the pre-equilibrium triton to 3 He 
ratio is also a sensitive probe to the density dependence 
of nuclear symmetry energy. 

Isospin effects on cluster production and isotopic ratios 
in heavy ion collisions have been previously studied using 
either the lattice gas model ^lj or a hybrid of IBUU and 
statistical fragmentation model These studies are, 

however, at lower energies than considered here, where 
effects due to multifragmentation as a result of possible 
gas-liquid phase transition may play an important role. 
Except deuterons, both tritons and 3 He are only about 
one per event in heavy ion collisions at 80 MeV/nucleon. 
The number of other clusters such as the alpha particle 
is not large either 0,0]. In this case, the coalescence 
model is expected to be a reasonable model for determin- 
ing the production of light clusters from heavy ion colli- 
sions. Furthermore, the effect obtained in present study 
will be enhanced if other clusters are emitted earlier as 
the isospin asymmetry of the residue is increased. 

In conclusion, using an isospin-dependent transport 
model together with a coalescence model for light clus- 
ter production, we have found that the nuclear symme- 
try energy affects significantly the production of light 
clusters in heavy-ion collisions induced by neutron-rich 
nuclei. More deuterons, tritons, and 3 He are produced 
with the stiff nuclear symmetry energy than the soft nu- 
clear symmetry energy. This effect is particularly large 
when these clusters have lower kinetic energies. Also, 
the isobaric ratio t/ 3 He, especially for higher energy tri- 
tons and 3 He, shows a strong sensitivity to the density 
dependence of nuclear symmetry energy. It is further 
found that light clusters production is not sensitive to 
the isospin-independent part of nuclear equation of state 
and the in-medium nucleon-nucleon cross sections. The 
study of light clusters production in heavy ion collisions 
induced by neutron-rich nuclei thus allows us to extract 
useful information about the density-dependence of nu- 
clear symmetry energy. 
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